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Abstract 



We reanalyse the cell theory of Lennard-Jones and Devonshire and find 
that in addition to the critical point originally reported for the 12 — 6 potential 
(and widely quoted in standard textbooks), the model exhibits a further crit- 
ical point. We show that the latter is actually a more appropriate candidate 
for liquid-gas criticality than the original critical point. 
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I. INTRODUCTION AND BACKGROUND 



Lennard- Jones-Devonshire (LJD) cell theory [0-0] is a lattice approximation to the liq- 
uid state. Historically it was the prototype microscopic model to predict the location of 
the liquid-vapour critical point of a simple fluid and, as such, is widely employed as an 
introduction to mean field theory in a number of standard texts and papers [[§-§] . It is still 
in use today for research applications 

Within the framework of the model, particles are considered to be localised in singly 
occupied "cells" , centred on the sites of a fully occupied lattice (of some prescribed symme- 
try), within which they move independently. For simplicity, cells are assumed to be spheres 
of volume v = V/N (the inverse number density) and radius s. A particle is considered to 
interact with its c nearest neighbours "smeared" around the surface of a further sphere of 
radius a concentric with the cell. The volume of this "interaction sphere" is related to the 
cell volume by 

a 3 = 1 v, (1.1) 

where 7 is a lattice-dependent constant, chosen so that for a primitive unit cell of volume v 
the lattice parameter will be the radius of the interaction shell (see figure [j](a)). 
The Helmholtz free energy of the model is given by JIJ 

NE n 

A = -Nk B T\nv f a c + . (1.2) 

Here E is the "ground state energy" - the energy per particle if all occupied their lattice 
sites; a c is a constant "communal entropy" term[], which attempts to account for the entropy 
lost due to the localisation of particles within cells; and Vf is the "free volume", 



I e -(E(r)~E )/k B T dv ^ (L3) 
Jv 



with E(r) the "cell potential" the interaction energy of a particle at a position r within its 
cell. 

Figure 0(b) shows the cell geometry. The radial symmetry of the system allows E(r) 
to be written in terms of E(r), where r is the radial coordinate of the particle's position 
within the cell. Simple trigonometry yields the separation R between the particle at point 
P within the cell and an element dA of the shell. If c neighbouring particles are assumed 
to be smeared over the interaction shell and the interparticle potential is given by u(R), the 
total energy of the particle E(r) is given by 

E(r) = c ^ U ^ dA , (1.4) 
Ana 2 

where, for a prescribed cell size, the density enters through the limits on R. One would 
normally perform this integration numerically with further numerical integration to calculate 
the free volume according to eq. |1.3| for a given choice of number density p = v~ x . 



1 Reference [|J assigns <r c = e, although the exact value chosen for a c has no impact upon phase 
coexistence so long as it is independent of volume. 
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L JD used the model described above to calculate the phase diagram for the 12-6 Lennard- 
Jones (LJ) potential, 



V 12 a 6 



u(r) = 4s — - - . (1.5) 



They assumed a face centred cubic lattice structure with coordination number c = 12 and 
geometrical constant 7 = y/2. The value of E(r) (eq. |L^) was calculated using the potential 



of eq. 1.5 , However for the calculation of the ground state energy Eq (appearing in eq. |P| 



they used the potential 



V 12 . a 6 



u(r) = As — - 1.2- . (1.6) 



Here the attractive part is increased by 20%, with respect to eq. [L^, which was apparently 
motivated as representing the effects of next and higher nearest neighbour interactions. 

For this system, LJD found a critical point at ksTcje ~ 1.3, P c /e ~ 0.6, p c ~ 0.56. These 
results are to be compared with the consensus arising from a number of simulation studies of 
criticality in the LJ fluid, the most recent and sophisticated of which |jTD[ quotes ksTc/e = 
1.3120(7), P c /e = 0.1279(6), p c = 0.316(1). That the cell theory apparently predicts the 
critical temperature of the LJ fluid to high accuracy has for many years been regarded as 
its major success ||. Indeed cell models remained the main method for calculation of liquid 
equations of state up until the advent of density functional theories [El. 



II. METHOD AND RESULTS 

We have obtained the cell theory phase diagram of the 12-6 LJ potential in order to 
compare with the results of LJD. Our study employed the same model parameters as used 
in the original study and was executed as follows. 

Using the Romberg method |J, numerical integration was carried out to calculate the 
cell potential E(r) and thence Vf and H (cf. eqns p..2| - |1.4| ). The equation of state in the 



space of the pressure P, volume per particle v, and the temperature T was obtained by 
using a golden section search |9j to minimise the Gibbs free energy per particle 

g(v) = A/N + Pv , (2.1) 

with respect to v across a spread of points in P, T space. Within this scheme, first order 
phase transitions were located by looking for jumps in v as P and T were variedF]. Critical 
points were estimated as the first point along a phase boundary where no distinct coexistence 
volumes could be found. 



2 Since the golden section search is vulnerable to finding local minima in g(v), we supplemented 
this approach with a Newton- Raphson root-finding algorithm |J, which solves for the coexistence 
conditions g{v\) = 5(^2), P(vi) = P{ v 2)- 
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The results of implementing this procedure, are shown in figure |^. We find a critical 
point at ksTcje = 1.354, P c /e = 0.53, p c = 0.57 which indeed appears to duplicate the 
original LJD critical point However, our full phase diagram shows two lines of first order 
transition, meeting at a triple point, both terminating in critical points. 

In attempting to identify the phases appearing in figure |^, it should be borne in mind 
that the lattice based character of the cell model precludes an accurate representation of 
non- crystalline phases. Accordingly it is not possible to assign unambiguously the physical 
nature of the phases. Nevertheless a plausible identification is suggestedbj the phase diagram 
itself which, in the region of the triple point, is reminiscent of that for a simple fluid such 
as Argon. Above the triple point temperature in such a system, an increase in pressure 
drives the system from a gas, to a liquid to a solid. That the same phase assignments are 
reasonable in the present model is supported by an analysis of their (radially symmetric) 
cell potential E(r), evaluated at the triple point (fcgTc/e = 1.02, P c /e = 0.025), which we 
now address. 

Figure 0(a) shows the cell potential E{r) expressed relative to the ground state energy 
Eq for the low density phase. This quantity is three orders of magnitude smaller than ksT 
across the whole cell, so the low density phase can quite reasonably be interpreted as gas-like. 
By contrast, the cell potential for the high density phase (figure 0(b)) exhibits a very steep 
minimum at r = 0, the centre of the cell. Away from the centre, the particles experience the 
hard core repulsion of their neighbours and the potential energy becomes large and positive. 
Thus the particle is strongly confined to its lattice site, and the phase can be considered 
solid-like. The cell potential energy for the intermediate density phase is shown in figure 
0(c). Its form is similar to that of the gas-like phase though the scale is different-there is 
a broad "hump" in the middle of the cell, of height around ksT. Thus the particle is still 
reasonably free within the cell and it seem reasonable to interpret this as a liquid-like phase. 
With these phase assignments 0, the critical point found by LJD is between a solid-like 
phase and a liquid-like phase, whilst the other, lower-temperature critical point terminates 
a liquid-gas transition at ksT c /e = 2.17, P c /e = 1.221, p c = 0.104. 



III. DISCUSSION AND CONCLUSIONS 

We have shown that the original LJD cell theory critical point for the 12-6 potential 
terminates a line of pseudo solid-liquid coexistence, rather than a liquid-gas transition 
as originally postulated. Moreover, we find that the model exhibits a further first order 
transition, which we believe is a more appropriate candidate for the liquid-vapour line than 
that suggested by LJD. It seems likely that LJD overlooked the transition we have found 
because at the time of their study, numerical integration was extremely labour intensive, 



3 Hill Q describes the cell theory without the modified ground state term used by LJD, but quotes 
the critical parameters corresponding to a modified system defined in ref. [12|. For the unmodified 
system we find cell potentials similar in form to those of fig. so we retain our previous phase 
designations. However the critical temperature of the "solid- liquid" -like critical point found by 
LJD is reduced to ksTcje = 0.82. The effect on the liquid-gas transition is much less pronounced. 
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rendering it difficult to attain the rather low densities at which the gas-like phase appears. 
Presumably also, the proximity of their critical point temperature to the "real" LJ critical 
point (estimated at the time from experimental data for Argon), made it tempting to make 
the identification they did. In hindsight this agreement can only be regarded as coincidental; 
one would expect a crude mean field theory such as the cell model to considerably overes- 
timate the critical temperature due to the neglect of fluctuations. The new critical point 
reported here does indeed do so. 

Of course a critical point between a solid-like and liquid-like states (as found by LJD) is 
an artifact of the cell structure of the model, which imposes long ranged orientational and 
translational order on the liquid where none exists in reality. The inability of cell theory to 
accurately represent disordered phases has long been appreciated and several attempts have 
been made to ameliorate this and other shortcomings of the approach. These include zero- 
and multiple-occupancy of cells ("hole" theories) 0,0; calculations of interactions with 
second- and higher nearest neighbours (see eg. ref. flT2"|); use of numerical integration to give 



more accurate cell shapes and free volumes [13]; and differing methods for calculating the 



cell potential (see eg. |/J). However, none of these extensions fully address the underlying 
lattice approximation which renders suspect any cell theory equation of state for a liquid - 
by definition, a disordered state without a lattice. Even hole theories will only be able to 
represent a liquid of appreciable density as a highly defective crystal. 

Finally we briefly point out implications of our findings for other recent work on cell 
theory. In ref. J?| a modified cell theory is introduced with the specific intention of engi- 
neering a freezing transition. This authors state that no such transition exists in the LJD 
cell theory, in the apparent belief that the critical point of the bare model is liquid-gas in 
character. 

In other recent work (which in fact motivated the present study), a cell model has been 
employed to study the phenomenon of the second critical point which is observed in pure 
fluids interacting via so-called core-softened potentials M. In such systems liquid-liquid 
phase separation occurs in addition to liquid-gas coexistence. Clearly in view of our finding, 
that even the basic cell theory for the 12-6 potential exhibits two critical points, the results 
of such studies should be interpreted with caution. We intend to report in depth on this 
matter in a future communication. 
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FIGURES 




FIG. 1. (a) A cell (dotted circle) of volume v and radius s is centred on a lattice site with 
primitive unit cell (dashed hexagon) of the same volume. The underlying lattice is shown with 
dot-dashed lines, with lattice sites as vertices. The nearest neighbours are assumed to be "smeared 
out" over the interaction shell (see text), of radius the lattice spacing a. (b) The geometry for the 
interaction of a particle at position P within the cell with an element of the shell dA. Given the 
radial coordinate r of the particle, the angular coordinate 6 of the element and the shell radius 
o, it is simple to calculate the separation R between particle and shell element; from this, the 
interaction can be integrated around the shell. 
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FIG. 2. (a) Projection of the cell theory LJ potential phase diagram in the (p,T) plane 



with p 



. Dots show density points found using the golden section search. Solid lines show 



coexistence densities, (b) (P, T) projection of the same phase diagram. Solid lines indicate lines of 
first-order phase transition, filled circles indicate critical points. Note that the lower temperature, 
higher pressure critical point is that found by LJD. 
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FIG. 3. Radially symmetric cell potential E(r) expressed relative to the ground state energy 
Eq, for the modified cell theory with LJ potential, calculated at the triple point from figure § 
(a) Potential for the low density phase (p/a 3 = 0.026). (b) Potential for the high density phase 
(p/a 3 = 0.869) (c) Potential for the intermediate density phase (p/cr 3 = 0.171), plotted from the 
centre of the cell, r/sa = 0, to the edge, r/sa = 1, where s is the radius of the cell at that density. 
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